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Abstract 

The acoustic response of a rigid frame porous slab with a periodic set of inclusions is cal- 
culated by use of a multipole method. The acoustic properties, in particular the absorption, of 
such a structure are then derived and studied. Numerical results together with a modal anal- 
ysis show that the addition of a periodic set of high-contrast inclusions leads to quasi-modes 
excitation of both the slab and the gratings, and to a large increase of the acoustic absorption 
of the initial slab, this being partly due to the quasi-modes excitation. 
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1 Introduction 

This work was initially motivated by the design problem connected with the determination of the 
optimal profile of a continuous and/or discontinuous spatial distribution of the material/geometric 
properties of porous materials for the absorption of sound. The equations that model the acoustic 
wave propagation in a macroscopically inhomogeneous rigid frame porous medium were derived in 
|10| . Acoustic properties of porous materials (foam) suffer from a lack of absorption particularly 
at low frequency. The usual way to solve this problem is by multi-layering |26| . 

In [21], the authors considered the reflection of a plane acoustic wave by a porous slab that 
presents a periodic set of pits. The medium is homogenized and its behavior is described in [S]. 
This leads to a drastic increase of the absorption coeficient at low frequency. In [27] the authors 
considered the transmission of an acoustic wave through a porous medium in which randomly- 
arranged metallic rods are imbedded, converted, by a procedure called ISA/3, into an equivalent 
homogeneous medium which exhibits decreased transmission and increased absorption. 

Periodic arrangements of either surface irregularities or volume heterogeneities usually lead to 
energy entrapment either at the surface or inside the structure, respectively, this being strongly 
linked to mode excitation, and to an increase of the absorption coefficient (first noticed by Wood 
|31j and partially explained by Cutler [7]). The particular properties of such structures have 
been studied in mechanics, with application to composite materials [5] 1171 i n optics, initially 
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motivated by the collection of solar energy [8l [18] , with applications to photonic crystals [HI [33] , 
in electromagnetics, with application to so-called left-handed materials [2S], in geophysics, for the 
study of the "city-site" effect [31 E] . The properties of such structures are now studied to create 
band-gaps for elastic or acoustic waves (phononic crystals [22, 23, 30J), but have only recently been 
used for the design of sound absorbing or porous materials [21 [2TJ [22] . 

Herein, we study the influence on the acoustic absorption of the introduction of a periodic set 
of fluid-like circular cylinders into a macroscopically-homogeneous porous slab (the porosity being 
homogenized in the equivalent fluid model). 



2 Formulation of the problem 

Both the incident plane acoustic wave and the slab are assumed to be invariant with respect to the 
Cartesian coordinate x%. A sagittal x\ — X2 plane view of the 2D scattering problem is given in Fig. 

m 




Figure 1: Sagittal plane representation of the configuration of plane wave solicitation of a d-periodic 
porous slab with fluid-like inclusions (of radius R) within a porous fluid-like slab. 

Before the introduction of the cylindrical inclusions, the slab is made of a porous material (e.g., 
a foam) which is modeled (by homogenization) as a (macroscopically-homogeneous) equivalent fluid 
M . Another equivalent fluid medium M 2 occupies each cylindrical inclusion. In the sagittal plane, 
the j'th cylinder is the circular disk Q 2 U) ■ The host medium M° occupying the two half spaces Q ± 
is air. Thus, we are dealing with a macroscopically-inhomogeneous slab, the heterogeneity being 
periodic in the x\ direction with period d. 

The upper and lower flat, mutually-parallel boundaries of the slab are T a and T^. The X2 
coordinates of these lines are a and b, the thickness h of the slab being h = a — b. The circular 
boundary of fi 2 co is ^(j) ■ The center of the j = disc is at the origin O of the laboratory system 
OX1X2X3. The union of f2 + and S7 - i s denoted by Oq. 

The wavevector k* of the incident plane wave lies in the sagittal plane and the angle of incidence 
is 9 l measured counterclockwise from the positive x\ axis. 
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3 Wave equations 

We designate a total pressure and wavenumber by the generic symbols p and k respectively, with 
p = p 0± , k = k° in Q ±,p = p 1 , k = k 1 in J7i and p = p 2 ^\ k = k 2 in Q 2 (j). 

Rather than solving directly for the pressure p(x, t) (with x = {x\, X2)), we prefer to deal with 
p(x,Lo), related to p(x, t) by the Fourier transform: 

/CO 
p(x,u)e~ lu;t du . (1) 
-00 

Henceforth, we drop the u in p(x, to) so as to designate the latter by p(x). This function satisfies 
the Helmholtz equations 

[A + (A; m ) 2 ] p(x) = ; x G fl m , m = ± , 1, 2 . (2) 
In an equivalent fluid medium [TU], the compressibility and density take the form 

1 7 P 



7 -(7-1) l + i^-G(PrV 



a; 



P?C^ r1 ^ (3) 



where w c = — — — is the Biot's frequency cut, 7 the specific heat ratio, Pq the atmospheric pressure, 

P/Ooo 

Pr the Prandtl number, pf the density of the fluid in the (interconnected) pores, <j) the porosity, 
the tortuosity, and a the flow resistivity. The correction functions G(Pr 2 cj) [2], F(u) [20] are 
given by 

G(Pr 2 w) = /^JpX, 

where 7/ is the viscosity of the fluid, A' the thermal characteristic length, and A the viscous char- 
acteristic length. 

The incident wave propagates in Q + and is expressed by 

p*(x) = A 1 exjp[i(k\xi - k\(x2 - a)) , (5) 

wherein k\ = —k°cos9' 1 , k\ = A;°sin#* and A 1 = A 1 (lo) is the signal spectrum. 

The new feature, with respect to the canonical case considered in [12J , is the transverse period- 
icity of Uj 6Z £} 2( ,-). 

Owing to the plane wave nature of the incident wave, and the periodic nature of Uj g z^2(j)> one 
can show that the field is quasi-periodic (Floquet theorem), i.e., 

p{( Xl + nd, x 2 )) = p({x 1 ,x 2 ))e [k i nd ; Vx G M 2 ; Vn G Z . (6) 

Consequently, it suffices to examine the field in the central cell of the slab which includes the disk 
O 2 (o) in order to obtain the fields, via the Floquet relation, in the other cells. Henceforth, we adopt 
the simplified notation: 0,2 : = J7 2 (o) , T := Fr \, p 2 = p 2 ^. 
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4 Boundary and radiation conditions 

Since M° and M 1 are fluid-like, the pressure and the normal velocity are continuous across the 
interfaces F a and IV 

p 0+ (x)-p 1 (x) = 0,Vx G r a , (7) 
(pT 1 9nP 0+ (x) - {p l y l d n p\ X ) = 0, Vx G r a , (8) 
A(x)- P 1 (x) = o, Vxer b , (9) 

(p )" 1 d n p ~ (x) - {p l )~ l c^ 1 (x) = 0, Vx G T b , (10) 

wherein n denotes the generic unit vector normal to a boundary and d n designates the operator 
d n = n • V. 

bmce M 1 and M 2 are fluid-like, the pressure and normal velocity are continuous across the 
interface T: 

p 2 (x)- P 1 (x) = o,Vxer, (li) 

(pT 1 d nP 2 (x) ~ (pT 1 d nP\*) = 0, Vx G r . (12) 

The uniqueness of the solution to the forward-scattering problem is assured by the radiation con- 
ditions : 

p 0+ (x) — p*(x) ~ outgoing waves ; |x| — > oo , X2 > a , (13) 
p° (x) ~ outgoing waves ; V |x| — > oo , %2 < b . (14) 

5 Field representations 

Separation of variables, the radiation conditions, and the Floquet theorem lead to the representa- 
tions: 

~° + (x) = [e^ (x2 - a) 5 p0 + R p e ik ^ (x2 - a) ] e lk ^ x \ Vx G , (15) 
p o- ( x ) = TpeK^i-^^-b)) j Vx G Qq _ ; (16) 
2pir 



p 



wherein S p o is the Kronecker symbol, k\ v = k\ H — — , k® p = \J (k ) 2 — (ki p ) 2 , with 5ft (k^p) > and 

The field in the central inclusion, with [r = (r, 6)}, takes the form 

p\v) =Y,CiJi (k 2 r) e iw , V(r,0) G ft 2 , (17) 

wherein J; is the Z-th order Bessel function. 

It is convenient to combine Cartesian coordinates (£1,0:2) and cylindrical coordinates (r,6) to 
write the field representation in Q\ . The latter takes the form of the sum (by use of the superposition 
principe) of the diffracted field by the inclusions p^(x) and of the diffracted field in the slab p*(x). 
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Because of the quasi-periodic aspect of the configuration, the diffracted field in the slab can be 
written in Cartesian coordinates as: 

PlW = E (fiv e ~ lklpX2 + ffp^ 2 ) e[klpXl • ( 18 ) 

pGZ 

We transform Cartesian to polar coordinates by means of: 

x\ = r cos (9) , k\p = k 1 cos (9 P ) , 
'4 P = kl sin {Op) > 



X2 = r sin (9) , k\ = k 1 sin I 



so that 



Use of the identity- 



leads to 



(19) 

J(r) = £ ( f-e^rco 8 (e+e p ) + ^yfcVco^-e^ (20) 
i>ez 

ifcirco S (0-e p ) = £ i m J m (kr)e lm{6 - 6 ^ , (21) 



e 



rfto = E E im (/r/ mep + /iV*) e 

m£Z p£Z 
m£Z pGZ 



„lm0 



(22) 



wherein J" p = i m e lm9 p and J+ p = i m e " lm ^. 

Let us now introduce: i) (r J (P), 9 J (P)), the polar coordinates of a point P in the system linked 
to the jth cylinder whose center is at the origin Oj , and ii) (rj , Oj ) the polar coordinates of the 
origin Oj in the polar coordinate system linked to the Zth cylinder. 

In the general case, the field diffracted by the inclusion appears as the sum of the fields diffracted 
by all the inclusions: 

= E E B3H ^ ( fcV ) eimeJ > ( 23 ) 

jez mez 

wherein Hm is the first-kind Hankel function of order m. 
Using Graf's formula [lj for the Hankel function leads to 

= E E E 4! } {k'r) e*""fz£. fl (*V ) j^M^ 
jezn{o} mez <jez 

+ ^ (jfci r ) e [me , for R<r<(d-R) . (24) 

meZ 

In the case of gratings, 7g = j X cZ, #q = it if j < or 9$ = if j > 0, and the quasiperiodicity 
implies that the multipole expansion coefficients relative to the jth cylinder of the grating are given 
by BL = e ik iB° = e ik iB m , Vm E Z. 
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The field expansion in the vicinity of the central cylinder then takes the following form (in 
agreement with 



Pd 



lei 



+ Y,Ji{k 1 r)e lW Y, S i-rnB m , for R<r<(d — R) . (25) 



lei 



raei 



with S l = Y {k l jd) [e lk ^ d + e~ ik ^ d . 

To derive an alternative form of (I25p in Cartesian coordinates, it is more convenient to start 
from Green's theorem 



Pd( r ) = 

wherein G is the Green's function 

G(r - r 

or, with [r s = (r s ,6 s )] and for R < x 2 



^(r s )G(r-r s )-^(r-r s )p(r s 



E 



2d £-t h 



1 

2p 



= ifci p (xi -x\ )+ifc2 P I X2 -X2 1 



fii _ s\ . _5_ _J_ \k\ p X\-V\k-i v x?. -\k\-px\-\kzpx\ 



l V -« 1 

pel 

J_ _J_ e ,iki p x 1 +ik 2 pX 2 e -ir s cos(6» s -6» p ) 
2d Z-L k 2p 



so that with the help of identity (I2ip . we are led to 

' E 1 

with 

r2?r 



Z = - / fcitafl ^ (-i) m J m (kR) e lm ^- e) Y 
Jo rnel lei 



B t Hf R + Y Si_ n B n ji(kR) 



r-2n 

+ / kRdB Y H)" 1 ^ [kR) e lm ^- e) Y 



1Z6» 



BtHf R + Y S l~n B nJl{kR) 



nei 



AW 



-ITT 



E<- 



i) m e lm H.R 



J m (kR)Hg\kR) - Hg\kR)j m {kR) 



Y(-ne [m9p B m , 



(26) 



(27) 



(28) 



(29) 



(30) 



wherein we have used the fact that J- m (kR)H$ (kR) - H$ (kR) j_ m (kR) 



2i 
irkR 



Proceeding in the same way for —R>X2 gives: 



pi(x) = — V J_ e ifcl « ,xl - ifc * xa x J , 



wherein 



o 



-m 



kRdd ("i) m ^m (kR) e lm ^ +e) Y 

m&Z I 

kRdO Y {-i) m j m {kR) e [m ^ +e ^ Y 



leZ 



m&Z 



in 

m — \raQ- 



lez 



BiHf R + £ Si- m B m ji(kR) 

m&Z 

BiHf R + £ Si- m B m Ji{kR) 



AW 



m£Z 



AW 



"kR 



J- m (kR)H$(kR) - H$(kR)j- m (kR) 



B r , 



^(-i) m e lmS »B„ 



The field diffracted by the inclusion is expressed in Cartesian coordinates by 

pez lez 

where the signs + and — correspond to X2 > R and x^ < R respectively, and 



K+ 

pm 



2( i) m AmOr, ts- _ ^( v „-im0 D 



imvp jy— 

dk\ v ' pm 



dk\ v 



with 9 P such that k l e Wp = k lp + ikL, [U [29]. 



(31) 



(32) 



(33) 



(34) 



6 Application of the continuity conditions 

Here we consider only the equations of continuity across the interfaces T a and IV The continuity 
conditions across T will be treated in section [71 



6.1 Continuity of the pressure field across T a 

From ([7]) we obtain 



2 p 0+ (( Xl ,a))e- ikllXl dxl - 1 2 p 1 {{x 1 ,a))e-' lkuXl dx 1 = 0, VI 6 Z . (35) 

d I d 

"2 J ~2 

Introducing the appropriate field representations therein and making use of the orthogonality rela- 
tion 



Mk in -k u )x ldxi = d5nu y^ n ) g 



gives rise to 



(36) 
(37) 



lez 
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6.2 Continuity of the normal component of the velocity across F 

From (TSl) we obtain 



(P°)" 1 ^ 2 p 0+ ((x 1 ,a))e- lfc «^ ( ixl 

d 

(p 1 ) _1 ^ 2 p 1 ((a;i,a))e- ifelia;i da;i = 0, VI 6 Z . (38) 



Introducing the appropriate field representation therein, and making use of the orthogonality rela- 
tion (f36|) . gives rise to 

- a° p 5 p0 + a° p R p + f^e^ - /+o£e ife ^ - £ #+i^e ifc 2> = , (39) 

wherein a„ = — Vz = 0, 1. 

p p l ' 

6.3 Continuity of the pressure field across T b 

From ([9]) we obtain 

2 p°~ (( Xl ,b))e- [kllXl dxl - / 2 P l {{ Xl ,b))e- {kllXl dxi = 0, VZ € Z . (40) 



d /_ d 

2 2 



Introducing the appropriate field representation therein, and making use of the orthogonality rela- 
tion ([36]) . gives rise to 

T P ~ /i"e- ife i 6 - /+e ife i 6 - £ ^Sje^ 6 = . (41) 

6.4 Continuity of the normal component of the velocity across 

From (|10p we obtain 

d 

ei 
2 

d 



(P°)" 1 ^ 2 p 0_ ((x 1 ,6))e- ifc » a;i ( ixl 

(p 1 ) -1 ^/^!, b))e- ikllXl d Xl = 0, VZ € Z . (42) 



Introducing the appropriate field representation therein, and making use of the orthogonality rela- 
tion (|36j) . gives rise to 

- a° p T p + /fXe- 1 ^ - f+ai^ + £ J^o^ = . (43) 



7 Determination of the unknowns 

Prom ([Ho]) . ([5U]) . (|41|) and ([4*5]) we get the expressions of / x ~ and in terms of B[. Introducing 
the latter into (TL8|) together with ([25]) leads to 



_ ^ p p ' (#+ e i(^i+fc^ 2 ) + ^- e i(fci P ^i-^ 2 ) 

pez zez ^* p 



«°) 2 - (« P ) 2 



pgz zez p 



x 



^+ e ifc 2p( a + fe ) e i( fc ip ;I 'i- fe 2p :I '2) _|_ j(- e -i k 2p( a + b ) e K kl P Xl+ i k 2p x 2)^j j (44) 



wherein 



-Dp = 2isin {k\ p L) \ {oP p ) 2 + (o^) J — 4Qp«pCos (k\ p L) , 

Flp " 1^ 6 * p0 ' ( 4 5) 

_ 2a p (a p - « p ) ifc i b 
t lp - e d p0 . 



'v 

To proceed further, we need to convert the cartesian form to the cylindrical harmonic form: 
p\r) = BiH\ l) {k\) e iw + £ J, (fcV) e iie £ ^_ mJ B m 

„iZ0 



+ E E ( J « + J i F i) 



lez pez 



+J2 J i ei/e E E (& n p - Bn > ( 46 ) 



wherein 



fa ) 2 - fa 1 ) 2 4f-iV-™ 
Qin P = V Pj ri n P x 4 cos (fc* p (a + b) + (Z + n) p ) , 

k 2 V U V d (47) 

flnp = A P n ~ x 1 j cos ((/ - n)9 p ) . 

Central to the multipole method are the local field expansions or multipole expansions around each 
inclusion [U QTJ [29] . Because p 1 (r) satisfies a Helmholtz equation inside and outside the cylinder 
of the unit cell, in the vicinity of the cylinder we can write 



p l (v) = Y,B l H i i 1) {k l r)e m + Y,Ji{k l r)e m A l . (48) 

zez 
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By identifying (|46j) with (jlSj). we find 

A = E ^— B ™ + E ("W + J i + E E (^mP - B ™ ■ (49) 



At this point, we account for the two equations (|12p and (|13f) . It is well-known that the coefficients 
of the scattered field and those of the locally incident field are linked by a matrix relation depending 
on the parameters of the cylinder only, i.e., 

R ^M&RWtfR) - 1 2 j l (k 2 R)J l (k 1 R) 

Bi = — : 7~j~\ ryn A i = R i A i i ( 50 ) 

7 2 J z (A; 2 i?)tf J (1) (fc 1 J R) - ^ x h\ \k 1 R)Ji(k 2 R) 

k j 

wherein 7 J = — r. Denoting B the infinite column matrix of components Bi, (1491) together with 
()50p may be written in the matrix form 

(I-RS-R(Q-P))B = RF . (51) 

with F the column matrix of mth element Jmp-^ip + JmpF\pi I the identity matrix, R the diag- 

onal matrix of component Ri and S, Q and P three square matrices of respective (m, q)th element 

Sm—qi ^ y Qmqp and ^ ' Pmqp- 
pgZ pgZ 

Remark: In case of a Neumann type boundary condition, the relation (]50p takes the form 

J* = ~3 klR) M := RfA t . (52) 
H\ 1] {k l R) 

8 Evaluation of the transmitted and reflected fields 

Once (jST]) is solved for Bi, V/ € Z, we can derive, from (f37|) and (jllj) . expressions for i? p and T p 
depending on Bi,\/l € Z, which, after introduction into (fT5|) . leads to the expression of the pressure 
field in O + : 

o +fxl = (*--)) ■ isin(^)((a M ) 2 -(a») 2 ) e i (^ + ^^- B )) 

P 1 J i ((a 0i ) 2 + (a H ) 2 ) sin (jfc«L) - ^"a * cos {kfL) 



peZ iez dfe 2p i (( Q 0)2 + (0,1)2) sin (fci p Lj - 2ai a o cos (A^i 

x a\ (-iaj sin (Z0 P + fc 2 » - aj cos (Z0 p + ^C^+^-a)) ; (53) 

and, after introduction into (fH)|) . leads to 
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p°-(x) 



-2a li a oi e i ( fc * a,1 ~ fc a i(a,a ~ 6) ) 



i((a 0i ) 2 + (a 

+ EE 

pGZ 2GZ 



sin 



2a H a 0i cos (k^L) 



elk 



x a 



2 P i ((a°) 2 + (a*) 2 ) sin [k\ p Lj - 1a\a® cos [k\ p L 
I (iag sin (Z0 P + fejpo) - a* cos (W p + fc^a)) eK*^-*^" 6 )) _ (54) 



9 Evaluation of pressure field in the slab 

We obtain from (|44jl the final expression of the field in Oi: 
2a 0i (ia 0i sin {k\ i {x 2 - b)) - a u cos {k\ i {x 2 - &))) 



^(x) 



i ((a 0i ) 2 + (a 1 *) 2 ) sin (jfe**L) - 2a u a 0i cos 



/e: 



P & i& dk *P i ((a£) 2 + («p) 2 ) sin ( fc 2 P L ) - 2a l a p cos ( fc 2p L ^ 
x cos{kl p x 2 + W p )e ik ^ L (a° p - a p ) 2 + ((a p ) 2 - (a°) 2 ) cos(A£ p (x 2 - (a + 6)) - W p ) . (55) 

Use of the continuity conditions on T leads to: 

„ N v-^ ^ 2 j l (k 2 R)j l (k l R)- 1 l j l (k l R)j l (k 2 R) , „ . lM 



Remark : The fields in ttj, j = + , 1, 0~ are the sum of i) the field in absence of the inclusions 
(whose expressions are the same as those in [12]) with ii) the field due to the presence of the 
inclusions. 

Remark : The field due to the presence of the inclusions, when compared with the Green's 
function as calculated in [12J in the case of a line source located in the slab, takes the form of the 
field radiated by induced periodic sources. The latter do not add energy to the system, but rather 
entail a redistribution of the energy in the frequency range of the solicitation. 



10 Modal analysis 

10.1 Modal analysis without inclusions 

In the absence of inclusions, the resolution of the problem reduces to: 



(q 11 - a 0i ) 



0i\ 



a 11 + a 0i ) 



-e ik ? L (a 1 



a ±l + a 
0i\ 



a 



T 
R 



ikfL ,1 



a 



0i\ 



(57) 
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The natural frequencies modes of the configuration are obtained by turning off the excitation. The 
resulting matrix equation possesses a non-trivial solution only if the determinant of the matrix 
vanishes, i.e. : 

i ((a 0i ) 2 + (a 1 *) 2 ) sin (k^L) - 2a u a 0i cos ($L) = . (58) 



1.01 r- 

1 v 
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Figure 2: Real part (left panel) and imaginary part (right panel) of the roots of the dispersion 
relation. 



Fig. [2] depicts the real and the imaginary parts of the solution c*{uo) = uj/{k\) of the previous 
equation with the mechanical characteristics used in the section [TT1 To solve the latter, we proceed 
as in |13j . The dispersion relation cannot vanish when the incident wave takes the form of a plane 
wave, because we must have 9? (&*) > k°, whereas k\ E [— ko, ko) in case of a incident plane wave. 



10.2 Modal analysis with inclusions 

For the slab with periodic inclusions, the problem reduces to the resolution of the linear system (|51|) . 
As previously, the natural frequencies of the modes of the configuration are obtained by turning off 
the excitation, embodied in the vector RF. The resulting equation possesses a non-trivial solution 
only if the determinant of the matrix vanishes: 

det(I-RS-R(Q-P)) = . (59) 

A procedure, called the partition method, for solving this equation, is not easy to apply because 
the off-diagonal elements of the matrix are not small compared to the diagonal elements. Even at 
low frequency (i.e. ?fc(k 1 )R « 1), Bj = Bf x (k l R) 2 + O ((^i?) 2 ) for j = {-1,0,1}, at least 
three terms should be taken into account. 

An iterative scheme can be employed to solve (|5ip and obtain an approximate dispersion rela- 
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tion. We re-write this equation in the form 



1-RilSo + Y, (Qup ~ p np) )) B i = 



Rl E Sl-m + E (Q lm 'P ~ Pimp) B m (1 - S m i) 



rneZ 



The iterative procedure for solving this linear set of equations is: 



p& 



i-Ri[s + Y, (Qiip - Pup) 



< B\ n) =R v 



i-rAso + Y, (Qup - Pup) 



E { J r P F i P + 4 F i P ) 



(61) 



+ E \ S l-m + E " J^mp) (1 " 5ml) 



, Vn € N* , 



(n) 

from which it becomes apparent that the solution of B} , to any order of approximation, is ex- 
pressed as a fraction, the denominator of which (not depending on the order of approximation), 
can become large for certain couples (ki p ,u) so as to make B^ n \ and (possibly) the field large for 
these values. 

When this happens, a natural mode of the configuration, comprising the inclusions and the 

(n) 

slab, is excited, this taking the form of a resonance with respect to B\ , i.e. with respect to a 
plane wave component of the field in the slab relative to the inclusions. As B^ is related to f^ p , 
f{~ p , T p and Rp, the structural resonance manifests itself for the same (ki p ,uj) as concerns the field 
in the slab and in the air. 



1 - Ri [ So + E (Qup ~ Pup) I = 1 " Ri [ E H o 1] 2 cos (k\jd) 

P ez J \j=i 



+ E Ahl 



pez dk zp Usm(kl p L 



ay 



+ K) 



2q°q;1 cos ( k\ p L 



(cos {k\ p (a + b) + 2W P ) (-1) 1 ((a° p ) 2 - (a*) 2 ) - e^ L (a° p - a*) 2 )) = . (62) 
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The latter equation is the sum of a term linked to the grating embodied in 1 — R[So with a term 
linked to the slab embodied in — R[ ]C P ez (Qiip ~ Flip)- This can be interpreted as a perturbation 
of the dispersion relation of the gratings by the presence of the slab. 



11 Numerical results 



'7-P0 



The ambiant and saturating fluid is the air medium (po = Pf = 1.213kg. m 3 , Co 

V p0 

with P = 1.01325 x 10 5 Pa and 7 = 1.4, 77 = 1.839 x 10 _5 %.m _1 .s _1 ). The infinite layer is 
1 x 10~ 2 m thick and filled with a polymer foam M 1 . The radius of the inclusions is constant 
equal to R = 2.5 x 10~ 3 m. We vary the center-to-center distance between each inclusion from 
d = 1 x 10~ 2 m to d = 2.5 x lCP 2 m. The inclusions are either filled with the air medium (which 
define the so-called type 2 sample) , with a melamin-foam (which define the so-called type 1 sample) 
or with such a material that the condition upon T is the Neumann one (which define the so-called 
type 3 sample). 

The medium M 1 is characterized by <p = 0.96, = 1.07, A = 273 x 10~ 6 m, A' = 672 x 10~ 6 m, 
(7 — 2843./V. s.tyi , while the nielairiine-foairi is characterized by <p — 0.99, c^oo = 

1.001, A = 

150 x 10~ 6 m, A' = 150 x 10~ 6 m, a = 12 x 10 3 N.s.m- 4 . The incident angle is 9 i = 0. 



The infinite sum ^ over the indices of the modal representation of the diffracted field by a 

M 

cylinder is truncated as ^ such that 



-M 



M = int (sR ^4.05 x {k x R) s + k x R^ + 10 . (63) 

On the other hand, the infinite sum ^ over the indices of the k\ v is found to depend on the 

pel 

frequency and on the period of the grating. We also use an empirical rule we have determined by 

p 

performing a large number of numerical experiments such that 

P =-p 

P = ^ + 50mt ^^CTTTC? j + 10 (_ ln ) ) + 5 . (64 ) 

In the latter equations int (a) represents the integer part of a. 

The developed form of the conservation of energy relation takes the form of 

1 = a + n + r , (65) 

with A, 1Z and T the absorption, hemispherical reflection and hemispherical transmission coeffi- 
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cients respectively, defined by 



2 p=—p 2 

flit) _ fSjt-t 



- p = —p ^ 



wherein p is such (hat ( k\ + — ( ^ + ^ > (A; ) 2 > ^jfej + . 
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Figure 3: Comparison of the absorption (left panel) and hemispherical transmission (right panel) 
coefficients between the cases of inclusions of type referenced by 1 (top panels), 2 (middle panels) 
or 3 (bottom panels) for various center-to center distances. 



The reflection coefficients are almost the same when compared with those as calculated in ab- 
sence of inclusions except for closed inclusions of type 3. The increase of the absorption coefficient 
is also due to the decrease of the transmission coefficient, as was previously noticed in [27]. In all 
the cases, the larger is the spatial period, the closer the coefficients are to those of the case of a 
homogeneous slab. 
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Figure 4: Comparison of the reflection coefficient with inclusions of type 3 for various center-to 
center distance. 



Fig. [5] compares the reflection, transmission and absorption coefficients for R = 2.5 x 10 



"m 



and R = 1 x 10 3 m type 3 inclusions when d = 10 x 10 3 m. The smaller the inclusions are, the 
closer the coefficients are to those of the case of the homogeneous slab. 
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Figure 5: Comparison of the reflection, transmission and absorption coefficients for R = 2.5 x 10 3 m 
and R = 1 x 10~ 3 m type 3 inclusions when d = 10 X 10~ 3 m. 



12 Discussion 

The addition of weak contrast inclusions (i.e. type 1 and type 2) leads to an absorption spectrum 
whose shape is close to that of a macroscopically- homogeneous slab. We also deal with an effective 
homogeneous porous slab, whose absorption coefficient is lower when the inclusions are filled with 
the air medium and larger when the inclusions are filled with an other, more absorptive, fluid-like 
porous medium. For fluid-like porous inclusions, the slope of the absorption coefficient is nearly the 
same as it is for the macroscopically-homogeneous slab, while for air inclusions (holes) the slope is 
larger. That is why at low frequency the absorption coefficient is lower than the one as calculated 
for the homogeneous slab and becomes larger for a higher frequency. 

The addition of high-contrast inclusions (i.e. type 3) leads to an absorption spectrum that is, 
on the average, larger than for the macroscopically-homogeneous slab, and which presents some 
additional peaks. The slope of the absorption spectrum is closer to that of the the macroscopically- 
homogeneous slab. Each peak leads to a supplementary increase of the absorption of the configu- 
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ration. The absorption can be multiplied by a factor of three at the location of some of these peaks 
when compared with the one of the macroscopically-homogeneous slab. Two types of peaks are 
distinguished. Each of them appears periodically for a specific period of the inclusions and can be 
associated with mode excitation of the whole configuration. 

Two peaks stand out for all center-to-center distances at the same frequencies near v = 65kHz, 
v = 133kHz. This corresponds to the excitation of quasi-mode 2 and 3 of the initial slab . A 
quasi-mode is a mode of the global configuration whose structure is close to the one of a mode of a 
sub-structure it is composed of. Quasi-mode 1 occurs for a frequency too low for it to be seen, and 
quasi-mode 3 is largely attenuated. The excitation of the quasi-mode of the slab is made possible 
by the spatial periodicity of the configuration. 

In the absence of the inclusions, the modes of a homogeneous slab filled with a fluid-like porous 
medium cannot be excited by an incident plane wave propagating in the air |15j . If cfet(E) = 
denotes a generic dispersion relation, then we can say that det(E) is close to zero (vanishes in 
absence of dissipation) only for an evanescent wave in the air medium, which cannot be excited by 
an incident plane bulk wave. 

The spatially-periodic configuration leads to a field representation, through the Floquet the- 
orem, that includes evanescent waves the air medium. Another explanation of the quasi-mode 
excitation relies on the fact that each inclusion acts as an induced cylindrical source. The re- 
sponse of an active cylindrical source radiating in the neighborhood of a homogeneous slab (see 
[12] ) enables a mode of this slab to be excited because some of the waves radiated by the source are 
evanescent and have the same structure (at resonance) as that of the evanescent wave associated 
with the mode. 

The other peaks appear in a periodic manner beyond the peak relative to the quasi-mode of 

the slab. Their periodicity is inversely proportional to the spatial periodicity d of the grating. This 

phenomena was already encountered in [13] and was attributed to a periodization of the quasi-mode 

2n 

due to the primitive reciprocal lattice vector ki = — ij.. This periodization is strongly associated 

a 

with the excitation of quasi Cutler- modes. 

The increase of the absorption is also due to an average (global) increase which does not take 
the form of additional peaks. This global increase can be explained by multipathing between each 
inclusion and/or by excitation of evanescent waves in the slab, thanks to the existence of the grating. 

Addition of high contrast inclusions also leads to an increase of the absorption of the slab, largely 
associated with a decrease of the transmission coefficient. This is mainly due to mode excitation 
of the whole configuration, whose structure consists of evanescent waves in the air medium (and 
in the slab), thus leading to an entrapment of the energy in the slab. The latter is dissipated by 
thermal and viscous effects. 

In all the cases, the closer the inclusions are, and/or the larger are their radii, the larger is the 
absorption. 

13 Conclusion 

We show that high-contrast, periodically- arranged, inclusions in a porous slab induce an increase of 
the absorption coefficient, mainly associated with a decrease of the transmission coefficient over a 
large frequency range, and in particular, at low frequencies (although higher than ~ 20kHz). This 
effect is due to mode excitation of the slab through excitation of the quasi-mode of the initial slab 
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(enabled) by the periodic inclusions, and to excitation of grating modes via multipathing between 
the inclusions. 

This increase of absorption is most noticeable for rigid frame porous inclusions in a large portion 
of the frequency range, and is less pronounced for air inclusions. 

The reflection coefficients are found to be of the same order as those in the absence of inclusions 
for low-contrast inclusions, and to be higher than those in the absence of inclusions for high contrast 
inclusions. The way of reducing the reflection is by acting on the surface geometry of the slab. A 
first approach can consist in the addition of a homogenized layer of double porosity [24\ . 
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A Numerical validation 

We validate the numerical implementation of the analytical calculation with the help of our Finite- 
Element code [16], in the case of a viscoacoustic medium (Fig. [6]). 

The medium filling f^o is the air medium (po = Pf, cq = * ). The infinite slab and the 

V Po 

inclusion are filled with viscoacoustic media whose relaxed characteristics are those of the porous 

medium when the dissipative aspects vanish, (i.e. p\ = — , ci = °° and no = — , 

P P ' V H 1.07 1-001 

C2 = - y== ) The dissipative aspects of M 1 and M 2 are described by constant quality factors 

over the frequency range of solicitation, such that Q 1 = 30 and Q 1 = 50. 

The numerical simulations are performed on a grid of 3000 x 400 nodes, whose grid spacing is 
equal to 1 x 10~ 4 m. The slab is 2 x 10~ 2 m thick and the center-to-center distance between two 
adjacent cylinders is 1.5 x 10" 2 m. The radius of each cylinder is R = 5 x 10~ 3 m. 

The signal spectrum of the solicitation is that of a Ricker pulse centered at vo = 100kHz. The 
incident angle is l = 0. 
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Figure 6: Comparison between the reflected and transmitted pressure field as calculated by the 
semi-analytical method and as calculated by the FE code. The top panels depict the reflected field, 
while the bottom panels depict the transmitted field for both a central location (left panels) and 
appendicular location (right panels) of the unit cells. 
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